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We analyze the influence of noise in transverse hexagonal patterns in nonlinear Kerr cavities. 
The near field fluctuations are determined by the neutrally stable Goldstone modes associated 
to translational invariance and by the weakly damped soft modes. However these modes do not 
contribute to the far field intensity fluctuations which are dominated by damped perturbations 
with the same wave vectors than the pattern. We find strong correlations between the intensity 
fluctuations of any arbitrary pair of wave vectors of the pattern. Correlation between pairs forming 
120° is larger than between pairs forming 180°, contrary to what a naive interpretation of emission 
in terms of twin photons would suggest. 
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I. INTRODUCTION 



The properties of the fluctuations and correlations in 
spatially extended patterns out of t her mo dynamical equi- 
librium were studied long ago in the context of hydrody- 
namic systems [^]J^ . More recently there has been a new 
surge of interest in the field of nonlinear and quantum op- 
tics since Lugiato and Castelli pointed out the existence 
of a purely quantum phenomenon in a spatial stationary 
dissipative structure ||] , the reduction of fluctuations be- 
low quantum limits in the difference between the inten- 
sities of the two Fourier modes of a stripe pattern. Since 
then, the properties of the fluctuations and correlations of 
stripe patterns in different nonlinear optical models have 
been widely studied [^-0]. Stripe patterns in these sys- 
tems appear as supercritical transitions, and close to the 
threshold for pattern formation the harmonics of the fun- 
damental wave vectors can be neglected. Therefore con- 
sidering the homogeneous mode and the two modes of the 
stripes for each field is enough and simplifies the problem 
allowing for an analytic treatment. The problem consid- 
ering the whole infinite set of transverse modes in optical 
stripe patterns has been addressed numerically in |^ . At 
difference with the few modes approximation, the con- 
tinuous problem makes evident the role of the so called 
soft modes in the fluctuations of the near field. However 
a very good agreement between the few modes approx- 
imation and the continuous treatment is still found for 
the correlations of the fluctuations of the far field mode 
intensities close to threshold. 

Due to the additional complexity, fluctuations and cor- 
relations in hexagonal optical patterns have been much 
less studied ||,|^. Furthermore, these two papers use 
an approximation in few modes to describe the pat- 
tern. However hexagonal patterns generally appear sub- 
critically, with a finite amplitude and therefore harmonics 
are not negligible even at threshold. Strictly speaking, an 
approximation in which only the homogeneous mode plus 
the six modes of the hexagons are considered is not fully 
justified. 



In this paper we will treat the case of a subcritical 
hexagonal pattern using a continuous model, i.e., avoid- 
ing any restriction to a reduced number of spatial modes. 
In particular, we consider an optical cavity filled with a 
nonlinear isotropic Kerr medium. This situation is de- 
scribed by the Lugiato-Lefever model |l^] in which a 
Turing instability was described for the self-focusing case. 
This pattern-forming instability leads to a hexagonal pat- 
tern in the 2D transverse plane. We find that the fluctua- 
tion of the near field are dominated by the two Goldstone 
modes associated to the translational invariance in the x 
and y directions, and by the soft modes arbitrarily close 
to them in infinitely large systems. We also show that the 
modes that dominate the near-field fluctuations do not 
contribute at all to the far-field intensity fluctuations. 
We identify the modes responsible for the far-field in- 
tensity correlations and we find: (a) strong correlations 
between arbitrary pairs of wave vectors of the pattern, 
but stronger between those forming a 120° angle; (b) 
anti-correlation between the zero wave vector of the spec- 
trum of fluctuations and any wave vector of the pattern. 
While the anti-correlation of the homogeneous mode with 
the off-axis wave vectors can be understood in terms of 
energy conservations, the common microscopic interpre- 
tation of the far-field intensity correlations in terms of 
emission of twin photons would naively suggest that the 
strongest correlation is between the wave vectors form- 
ing a 180° angle. In fact, the total transverse momentum 
conservation involve always at least four modes simulta- 
neously and give some hints about how the correlations 
should be, but does not identifies the pairs with stronger 
correlations. Our results here are obtained within a semi- 
classical approach in which specific features of quantum 
statistics are neglected. 

The paper is organized as follows: In Sec. El we de- 
scribe the model we are considering. In Sec. [II we lin- 
earize around the hexagonal pattern and describe the lin- 
ear response of the system to noise perturbations. In Sec. 
IV we discuss in detail the field fluctuations, and in Sec. 
^we describe the correlations of the field Fourier com- 
ponents. Finally, in Sec. ^ we give some concluding 



1 



remarks. 



II. DESCRIPTION OF THE MODEL 

The dynamics of the electric field inside an optical cav- 
ity with a self-focusing Kerr medium can be described, 
in the mean-field approximation, by a equation for the 
scaled slowly varying amplitude of the field E{x) IlOlJlT 



dtE 



^i2\E\^E^^{x,t), 



(1) 



where Eq is the input field, 6 is the cavity detuning, 
is the transverse Laplacian and ^ (x, t) is a complex Gaus- 
sian white noise with zero mean and correlations 

{^{x,t)C{x'.t')) = eS{x-x')S{t - t') 
(^(f,t)e(x^O)=0. (2) 

In the absence of noise Eq. (|]) has a homogeneous sta- 
tionary solution Eg given by 



^0 = ^.[1-^(2/.-^)], 



(3) 



where Is = \Es\^ . It is well known that the homogeneous 
solution (||) shows bistability for 6 > ^/S. We will restrict 
ourselves to the non-bistable regime < V^- 

A linear stability analysis of the homogeneous solution 
with respect to spatially periodic perturbations yields to 
the dispersion relation 



X{k) = -l± y^-{0^k^-6Is){0^k^-2Is), (4) 

where X{k) is the linear growth rate of a perturbation 
with wave vector k and k = \k\. The instability thresh- 
old is located at 1^ = 1/2 and the critical wave number 
is kc = \J —B -|- 2. For pump intensities above threshold 
the maximum linear growth rate is for wave vectors with 
modulus 



(5) 



At threshold ku = kc- Starting from the homogeneous so- 
lution and changing the pump intensity to a value above, 
but close to, threshold a hexagonal pattern with a wave 
number k close to ku arises (which as follows from Eqs. 
(H) and (^) depends on the pump intensity). The tran- 
sition is subcritical and the hexagonal pattern, once it is 
formed, is stable for values of the pump intensity within a 
quite large range which includes values below threshold. 
Typically, the hexagons appear oscillating even at thresh- 
old. The amplitude of the oscillations decreases decreas- 
ing the pump intensity until they become, in all cases, 
stationary. We are interested here in the properties of 
the fiuctuations and correlations of the stationary hexag- 
onal patterns. The oscillatory behavior of the hexagons 
is investigated in Ref. [O]. Due to its subcriticality, the 



hexagonal pattern has always a finite amplitude, even at 
threshold. The harmonics of the six fundamental wave 
vectors have also a significant amplitude so that they 
have to be included in the calculations to obtained re- 
alistic quantitative results. This is particularly relevant 
here due to the self-focusing effect which leads to a pat- 
tern with high peaks and a strongly anharmonic far-field 
(See fig. 0). 
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FIG. 1. (a) Near field intensity \E{x,y)\'^ and (b) power 
spectrum (far-field) \E(k)\'^ of a stationary hexagonal solu- 
tion. Note the presence in the far-field of high order harmon- 
ics, (c) Shows a cross section along a x axis of the intensity 
pattern and (d) a cross section along the ky axis of the power 
spectrum. 



III. LINEARIZATION AROUND THE 
HEXAGONAL PATTERN 



The stationary hexagonal pattern can be written in the 
form 



N 



ik^ (x—xo) 



(6) 



n=0 



where are complex coefficients, xq determines the 
global position of the pattern (we take = in the 
following), is the homogeneous mode and k^ for 
n = 1, . . . , A/" are the off- axis wave vectors of the hexago- 
nal pattern. Here we take TV = 90 which corresponds 
to considering up to the fifth order harmonics in the 
far field. The six fundamental harmonics have modu- 
lus kw For simplicity we first consider the case without 
noise. Linearizing Eq. (|]) around the stationary solution 
(^ we obtain the following equation for the fluctuations 
SE{x,t) =E{x,t) -Eh{x) 



2 



dt5E = -(1 + ie)5E + iV^6E 

+i2[2\EHfSE + EhEhSE*]. 



(7) 



As Eq. (Q) is a linear differential equation with periodic 
coefficients, a general bounded solution can be found un- 
der a Floquet form [|l3|: 



(8) 



where A{q^ t) are functions with the same spatial pe- 
riodicity than the stationary pattern Eh^ and therefore 
can be written as 



N 



(9) 



From (|^), (H) and we obtain for each perturbation 
wave vector q a set of linear differential equations for the 
time evolution of the Fourier coefficients dan{q^t). 



dtSan{q,t) 



= [-{l^ie)-i\k'^^q]^]SaMt) 

N N 

n — l-\-m {q,t) + 



1^0 m=0 
N N 



aiam[Sa-n^i^m{-q,t)r} (10) 



^=0 m=0 



where San-i^m{q,t) = Saj{q,t) with kj=k^- + k^. 
Introducing T^{q^t) = {Re[Sao{q^t)]^ Im[Sao{q^t)]^ . . ., 
Im[5aN{q^ t)], Re[5ao{—q, t)], . . ., Im[5aN{—q, Eqs. 
(|lO|) can be written as 



dt^{q,t)=M{Eh,q)^{q,t). 



(11) 



S(^, t) includes perturbations with -\-q and —q since they 
are coupled in Eq. (p!o|). The important point is that per- 
turbations with different q vectors are uncoupled. Any 
perturbation with a vector q' outside the first Brillouin 
zone of the hexagonal lattice defined by the wave vectors 
of the pattern k^ (See fig. ^ is equivalent to another one 
with a vector q = q^-\-k^ inside. 




FIG. 2. First Brillouin zone (dashed hexagon) of the hexag- 
onal lattice defined by the wave vectors of the pattern in the 
Fourier space. 

Therefore the values of q to be considered are only 
those inside half of the first Brillouin zone. In this way 
one finds a set of 4 (A^ + 1) eigenvalues and eigenvectors 
for each vector q except for the cases q = and q = k^/2 
where S(g, t) has only 2(A'+ 1) components. The eigen- 
values may be either real or complex conjugates and de- 
termine the linear response of the system when pertur- 
bations with k^^q wave vectors are applied. 

As the system is translationally invariant, Eh{x -\- xq) 
is also a stationary solution for any fixed xq, the two 
modes dxE^ and dyE^ are eigenvectors of M{Eh^q = 0) 
with zero eigenvalue. These neutrally stable modes of the 
linearized dynamics are the so-called Goldstone modes 
[|l|,0 and correspond to homogeneous perturbations that 
displace rigidly the pattern along the x or y directions 

dM^)=i^^nkLe'^-^ (12) 

n 

Eh{x, y) + xod^Eh = ^ o„(l + ikl^xo)e'^'^^ 

n 

^ a„e*^"(-+--f) = Eh{x + xo, y). (13) 



For the hexagonal pattern considered here the Goldstone 
modes have the profile shown in Figs. || and ^. In the 
near field they have a larger amplitude at the borders 
of the hexagonal peaks, where the gradient is larger and 
show a sharp transition from negative to positive values 
just at the center of the peak (Fig. ^). 




FIG. 3. Real part (left) and power spectrum (right) of the 
Goldstone modes dxEh (top) and dyEh. 
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FIG. 4. Close up of the real part of a peak of the hexagonal 
pattern (left) and the real part of the Goldstone mode dx^h 
(right). 

In a situation where the hexagonal pattern is stable 
the Goldstone modes are the only neutral modes while 
all the other eigenvectors have eigenvalues with strictly 
negative real part. 

At this point we can study the linear response to noise 
perturbations. Decomposing the noise term of (|l|) as it 
has been done for the field [Eqs. (||) and (^)] we obtain 
an extra term in the rhs of equation (O) 



dtt% t) = M{Eh, q)t{q, t) + S(g, t), 



(14) 



where ^{q,t) = [ReUk^^ + q,t% ImUk° + q,t% .., 

ImUk%+q,t% ReUk^-lt)],..., ImUk% - q,t)]Z- 
Then, we can write an Ornstein-Uhlenbeck process Q] 
for the amplitude Oi(<f, t) of each eigenvector of M{Eh^ q) 



(15) 



where Xi{q) is the i-th eigenvalue of M{Eh,q) (ordered 
according to the value of its real part, Re[Xi{q}] > 
Re[X,+^{q)]). vMt) = EtZVCr/m^m) is the 
noise expressed in the eigenvectors basis and C{q) is the 
matrix for the change of basis as obtained diagonalizing 
M{Eh^q). The coefficients 6an{q,t) are related to the 
amplitudes of the eigenmodes Oi{q^t) by 

Re[San{q,t)] = ^iei{q,t)C2n-ii{q) 
Im[San{q,t)] = ^iQiiq,t)C2ni{q) 
Re[dan{-q,t)] = T^iQ i{q, t)C2n -1+2 Ni{q) 
Im[6an{-q, t)] = T^iQi{q, t)C2n+2Ni{q)' (16) 

The noises Tjii^q^ t) are Gaussian and white in time and 
have cross-correlations 

{rii{q, t)ri*{q', t')) = '-D^jimt - mq- q') (17) 



{r]^{q.t)mW ^t')) = ^l^D,^M8{t-t')3{q-q'). (18) 



where 



4Ar+3 



k=0 
4Ar+3 



(19) 



The solution of the stochastic process is then 



e-^^^^>r]Ms)ds. (20) 



The average value of the eigenmodes amplitude is 
{Oi{q,t)) = 0. The correlations between the amplitudes 
of the eigenvectors are given by: 



*(9'.0) = '-£Mm_^ 

-4(A.(9-)+A*(9')) 



(1 _ e(^-(9)+^.*(9*'))*)5(g'_ q') (21) 



{eMt)ej{q',t')) 



eDijiq) 



-4(Ai(g-) + A,- (</')) 
(1 _ e(^*(9)+^i(9''))t)5(g _ q') (22) 



In particular, the time evolution of the mean value of 
the squared amplitude of the eigenvectors with non zero 
eigenvalue is 



For times much longer than a characteristic time ri{q) ~ 
— l/Re[Xi{q)] this mean squared amplitude reaches a sta- 
tionary value 



(mm = 



2- _ eDu{q) 



-8Re[Xi{q)y 



(24) 



Eq. ( |24| ) does not apply to the Goldstone modes as they 
have zero eigenvalue {Xo(q = 0) =0). Its time evolution 
is given by (from Eq. (pQ)) 



r 

Jo 



r]i{q,s)ds. 



(25) 



This is a purely diffusive motion which never reaches a 
stationary state. Its mean squared amplitude grows lin- 
early in time 



|eo(g = o,t)|2) = -Doo(9 = o)t. 



(26) 



/c=0 



Therefore, the linearization fails for times t ~ 1/e, when 
the amplitude of the Goldstone modes |0o(^ = 0)p 
reaches values comparable to one and nonlinear terms 
will come into play. 



4 



IV. FLUCTUATIONS IN THE NEAR FIELD 

Starting from a stationary hexagonal solution of the 
system equations without noise, a typical evolution of the 
pattern fluctuations when the noise is switched on, ob- 
tained from numerical integration of the nonlinear equa- 
tions Eq. (H) jl^], is shown in Fig. ^ 
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FIG. 5. Near field (left: Re[5E{x)], 5E{x) = E{x)-{E{x))) 
and far field (right: |^^(^)|, 5E{k) = E{k) - {E{k))) pattern 
fluctuations after switching on the noise. From top to bottom: 
t = 2, t = 200, and t = 2000. We have considered e = 10"^ 



From Eq. (23) we get that for short time all the modes 
Qi{q^t) have a similar mean squared amplitude which 
is proportional to the noise intensity and grows linearly 
with time 



(27) 



All the modes will contribute with a similar weight to 
the field fluctuations and therefore there is a complete 
lack of structure in the field fluctuations at short times 
as shown in Fig. jH] (top). 

As time goes on, following (p3[), the mean squared am- 
plitude of the modes does not grow linearly any more. 
It reaches a steady state value given by ( p^ which is 
larger for the modes which have a smaller decaying rate 



Re[Xi{q)]. According to the way the eigenvalues has 
been ordered, for a given the smaller decaying rate 
is Re[Xo{q)]. The eigenmodes with eigenvalue Xo{q) are 
the so-called soft modes |ll,[7|,|l6[| which are connected with 
the Goldstone modes and for which Ao(<f ^ 0) ~ —\q\'^ 
(see Fig. 1^). 




-0.04 E 

0,0 0.1 0.2 0.3 0.4 0.5 

FIG. 6. Left: Re[Xo{q)], darker color indicates smaller val- 
ues. The white line shows the first Brillouin zone. The 
center of the figure corresponds to the Goldstone modes 
{Re[Xo{q = 0)] = 0). Right: transverse cut of Re[Xo{q)] along 
qx (solid line) and qy (dotted line) axis. 

While the Goldstone modes correspond to neutrally 
stable homogeneous perturbations, the soft modes corre- 
spond to weakly damped long- wavelength perturbations. 
For systems with a finite size L the less damped of the soft 
modes are the ones with the smallest q vector allowed by 
the size of the system, \q\ = 27r/L. These modes have a 
decay rate proportional to and therefore a station- 
ary mean squared amplitude proportional to eL^. The 
shape of one soft mode is illustrated in Fig. |^. 



• • • • 



FIG. 7. Soft mode q 

Re[X{q^ {0,27V /Ly))]^- 




= (0,27r/L^ 
-0.0052. 



whose decay rate is 



Its profile is similar to a Goldstone mode but with a 
long- wavelength modulation on top of it. The wavevector 
of the modulation is precisely the wavevector q that iden- 
tifies the soft mode. While the Goldstone modes describe 
an overall rigid motion of the pattern, the soft modes add 
some distortion so the pattern moves in slightly different 
form at different locations. 

For the numerical simulations shown in Fig. |^, at 
t = 200 the mean squared amplitude of most of the modes 
has already saturated at small values, while the ampli- 
tude of the soft modes is, in average, just reaching the 
stationary value. Therefore, as shown in the figure (cen- 
ter) the field fluctuations are dominated by the soft and 
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Goldstone modes and short range spatial structures start 
to appear. The typical correlation length of the fluctua- 
tions is determined by the wavelength of the soft modes. 

According to ( p^ ) the mean square amplitude of the 
Goldstone mode keeps growing linearly in time and does 
not saturate in the linear approximation. The Goldstone 
mode amplitude only saturates due to nonlinearities at 
times of the order t ~ 1/e where the linear theory for 
fluctuations described in the previous section fails (this 
time is around 10^ for the simulation shown in Fig. ||) As- 
suming linearization is valid, at times t ^ L^, the mean 
squared amplitudes of the Goldstone modes (~ et, (|2^)) 
become much larger than the amplitude of any of the soft 
modes and dominate the fluctuations. In our system at 
time t ~ 2000 the profile of the field fluctuations is al- 
ready determined only by the Goldstone modes as shown 
in Fig. H (bottom). In the far field the largest fluctua- 
tions are those of the wave vectors of the pattern, while 
in the near field fluctuations show a long range structure 
with strong correlation over all the system size. Inte- 
grating for longer times will not change the profile of the 
field fluctuations, the amplitude of the Goldstone mode 
will eventually saturate but nevertheless will completely 
determine the field fluctuations. 

For systems with a larger system size, the soft modes 
will take a longer time to reach the stationary amplitude 
(and in fact its amplitude will be larger). However, pro- 
vided that < 1/e, the evolution will be basically the 
same, showing the three stages discussed before. If the 
system size is larger than then the stationary am- 

plitude of the soft modes and the Goldstone modes will 
be determined by the system nonlinearities and if L is 
large enough, both amplitudes may have similar sizes, in 
which case Goldstone and soft modes will contribute to 
the near field fluctuations even at t ^ oo and therefore 
the pattern of fluctuations may show a finite correlation 
length in the stationary regime. 



V. FLUCTUATIONS AND CORRELATIONS IN 
THE FAR FIELD 

In the following three subsections, we address the fluc- 
tuations in Fourier space of the field, intensity and mo- 
mentum. 



A. Field fluctuations 

In the Fourier space the field fluctuations are also dom- 
inated by the Goldstone and soft modes. From Eqs. ([l^ ) 
and (|l|), we can see how the Goldstone modes, which 
have the same wave vectors than the hexagonal pattern, 
induce opposite and very large phase fluctuations in op- 
posite Fourier components, which correspond to the rigid 
translation of the hexagons in the near field. The homo- 
geneous component of the field is not affected by the 



Goldstone modes. The soft modes do not have exactly 
the same Fourier components than the hexagonal pattern 
but they are very close, therefore their main contribution 
is to broaden the spots of the far-field fluctuations (see 
the far-field for the intermediate time in Fig. ^). 

B. Intensity fluctuations and correlations 

The intensity fluctuations of the far-field peaks are 

5I(k)=I(k)-h(k), (28) 

where I{k) = \E{k)\^ and Ih{k) = \Eh{k)\^. The corre- 
lation function of the far-field intensity fluctuations of a 
fundamental wave vector of the pattern, for instance ^3, 
with the far-field intensity fluctuations of any other wave 
vector k is given by 

C.(fcl,fc) = ^J£ai£i^. (29) 
^{\SI{kl)\^){\SIikW) 

From the numerical integration of Eq. (|l|),and averag- 
ing over two hundred realizations of the noise [|l5[ , we find 
strong correlations between the intensity fluctuations of 
all the modes of the pattern, not only among the fun- 
damental harmonics but also with the higher order ones 
(Fig. I). 




FIG. 8. Correlation function of the intensity fluctuations 
Ci{ks,k). The brightest spot takes the value 1 correspond- 
ing to the autocorrelation of ^3. Note the strong correlation 
with the other six fundamental wave vectors of the pattern 
and with the first rings of harmonics. The correlation decay 
as one consider higher and higher harmonics. Note also the 
ant i- correlation with the homogeneous as a black spot at the 
center. The grey background is around 

For the fundamental harmonics, the correlation 
Ci{klk^) = {SI{k^)6I{k^)), with z = 1,6, is larger for 
the modes forming an angle = ^ and smaller for ^ = |^ 
(Fig. [1^). Using a 6 modes approximation ||,|^ one finds 
that Ni + Ni^i - Ni^s - Ni^^ (being Ni the number of 
photons of the mode i ) is a conserved quantity for the 
interaction Hamiltonian related to momentum conserva- 
tion. This means that one should expect strong correla- 
tions among these sets of four Fourier modes. However 
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this reasoning can not reveal which modes within these 
sets are more correlated in pairs. So, the stronger cor- 
relation between those modes forming; a (9 = ^ angle, 
despite fulfill momentum conservation, can not be com- 
pletely understood in this terms. We also find strong 
anti-correlations between the intensity fiuctuations of the 
modes of the pattern and the homogeneous mode which 
are related to energy conservation (Fig. ||). 

Neglecting terms of order e^, the far-field intensity fiuc- 
tuations can be approximated as 



'i{C2m-liC2nj + C2n-liC2mk)) (32) 

{5E{kl)5E*{kl:)) = 
E, 34(i7Tip (1 - e(^^+^P*)A,(g = 0) 

{C2n-liC2m-lj + C2niC2mj + 

'i{C2m-liC2nj — C2n-liC2mk)) (33) 



SI{k) ^ 2Re[El{k)6E{k)], 



(30) 



where 5E{k) = E(k) - Eh(k). As Eh(k) = 
an(27r)^(5(^5J — k) we have to consider only pertur- 
bations that have the same wave vectors of the pattern 
{q = 0). Therefore, the soft modes {q ^ 0) do not con- 
tribute to the fiuctuations of the far-field intensity peaks. 
The Goldstone modes are indeed q = perturbations, 
but the intensity fiuctuations associated to the Gold- 
stone modes are SI{k) = 2Re[Es{kY dxEgik)] which from 
(^ and ( p!2[ ) is exactly zero. So, neither the Goldstone 
modes do contribute to the fiuctuations of the far-field 
intensity peaks. They only contribute to phase fiuctu- 
ations. Therefore, fiuctuations of the far-field intensity 
peaks have to be described by the other eigenvectors of 
M{Eh, q = 0). These eigenvectors have the same Fourier 
components than the hexagonal pattern and their eigen- 
values Ai(^= 0), (z = 1, . . . A^) are shown in Fig. 0. 




FIG. 9. Eigenvalues of the eigenvectors with q = 0. There 
is two zero eigenvalues corresponding to the Goldstone modes. 
Note also the symmetry of the spectrum with respect to the 
axis i?e[A] = — 1. 

Using Eqs. (^6|), ( |2l| ) and ( p^ ) we can compute the 
correlation function Ci{k^^k^) analytically 



Ci{k^, k^ 



where 



2Re[EUk'jEUk°J{SE{k°J5E{kl^)) + 

EH{Z)El{kl){5E{Z)5E*{kl)), (31) 



Fig. |10| shows the stationary value (t ^ oo) of ( |3l[ ) for 
the fundamental wave vectors of the pattern as a function 
of the angle 9 between them. The results obtained from 
the linearized theory are in very good agreement with the 
correlations obtained from the numerical simulations of 
the nonlinear equation (||). 
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FIG. 10. Correlations between the intensity fluctuations of 
the fundamental wave vectors obtained from the numerical in- 
tegration of equation (|^) (rhombi) , and analytically form Eq. 
(El) (dots). 



From Fig. y we can see that the most important eigen- 
vectors are those associated to the complex conjugate 
eigenvalues with Re[\2{q = 0)] = Re[Xs{q = 0)] = —0.25. 
This pair of eigenvectors give the strong correlation be- 
tween all the Fourier components and a strong anti- 
correlation with the homogeneous field. The excitation 
of this eigenvectors implies also anti-correlations between 
the homogeneous mode and the off-axis Fourier compo- 
nents of the pattern (Fig. O). 
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FIG. 11. Intensity fluctuations due to the eigenmodes with 
Re[X] = -0.25 



C. Transverse momentum fluctuations 

We finally address the fluctuations of the transverse 
momentum. Without noise the total transverse momen- 
tum of the pattern is P = \an\'^k^ = 0. The noise 
induce fluctuations given by SP = ^^2Re[a^San]k^' 
One finds that all the eigenvectors that contribute to the 
far-field intensity fluctuations strictly fulfill momentum 
conservation except for two modes with X2n{q = 0) = 
X2N-i{q = 0) = —2, which are exactly those symmetric 
to the Goldstone modes respect to the line Re[X\ = — 1 
(Fig. 1^). Therefore momentum fluctuations are deter- 
mined by the two modes with maximum damping which 
are shown in Fig. |l^, the on the top figure breaks the 
conservation of the Py momentum component while the 
one in the bottom breaks the conservation of Px. 



Some eigenvectors with Re[Xn{q = 0)] = — 1 are finally 
the responsible for the differences between the correla- 
tions in the fluctuations of the six firth harmonics. The 
typical profile of one of these eigenmodes is shown in Fig. 
|l^. The symmetry of the three peaks located at 120° 
yields to larger correlations between the fluctuations of 
the corresponding far field intensity peaks. 
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FIG. 12. Intensity fluctuations due to one of the eigen- 
modes with i?e[A] = — 1 



We note that the modes that determine the correla- 
tions of the far-field intensity peaks reach stationary val- 
ues in a much shorter time {t ~ —l/Re[X2{q = 0)] = 4) 
than the Goldstone and soft modes contributing to the 
near field fluctuations {t ~ 1/e ~ 10^ as discussed in 
the previous section). As these correlations are deter- 
mined by ^ = perturbations, these calculations can be 
performed in systems with relatively small size. Increas- 
ing the system size, may change the near field profile of 
the fluctuations, as discussed before, but nevertheless the 
dynamical evolution of the modes contributing to the far- 
field intensity fluctuations and their mean squared sta- 
tionary amplitude will be basically the same. 
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FIG. 13. Intensity fluctuations due to the two modes 
Re[X2N{q = 0)] = Re[X2N-i{q = 0)] = -2. 

Alternatively it is possible to see from the classical field 
theory that the momentum fluctuations are damped with 
a coefficient -2. A Lagrange density can be defined for 
Eq. (0) 



L = e2*[a|VE|2 + -{EE* - EE*) - 
i{EoE* - E^E) + 0\E\^ - \E\ 

Substituting 



(34) 



in (34) and integrating over all the space 
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we get a Lagrangian for the amphtude of the raodes 
and for the position of the rohs in the near field xq 

dna; - i2\an\^knXo) + i{Eoa^ - E^ao) + 

En l^^l ~ En En' Em ^^^^'^m^n+n'+ml • 

(35) 

xq is a cyclic coordinate, so its conjugate momentum 
dL 



dxo 



(36) 



is constant. Therefore P = — 2P. P is identically zero 
without noise. When noise is present momentum fluctu- 
ations should satisfy 



6P = -2SP^x{t). 



(37) 



where x{t) is a Gaussian white noise. Therefore, mo- 
mentum fluctuations have the maximum damping. The 
same damping coefficient for the momentum ffuctuations 
has been found by Gatti and Mancini from a few modes 
quantum formulation [51. 



VI. CONCLUSIONS 

We have analyzed the ffuctuations and correlations in 
a hexagonal pattern of a prototypical model in nonlinear 
optics. In the near ffeld ffuctuations are dominated by 
the neutrally stable Goldstone modes associated to the 
translational invariance as well as by the soft modes con- 
nected with them. The soft modes destroy the long range 
correlation in the ffuctuations, however in small systems 
these modes reach an stationary amplitude much earlier 
(and at a smaller value) than the Goldstone mode so that 
they are important only at intermediate times. At long 
times the ffuctuations are dominated by the Goldstone 
modes which correspond to rigid displacements of the 
overall pattern. For very large systems, both the Gold- 
stone and soft modes may have similar amplitudes and 
contribute to the ffuctuations. 

In the far-ffeld, the most relevant effect of noise are the 
intensity ffuctuations of the Fourier modes of the hexag- 
onal pattern. At ffrst order in noise intensify, these ffuc- 
tuations are not affected neither by the Goldstone modes 
nor by the soft modes. They are dominated by damped 
modes, so they reach stationary values in relatively short 
times. Their main characteristics are: i) strong correla- 
tions between the intensity ffuctuations of any arbitrary 
pair of the six fundamental wave vectors of the pattern, 
and also with their higher harmonics, ii) larger correla- 
tion between intensity ffuctuations of the Fourier modes 
forming 120° angles than between modes forming 180°, 
and iii) strong anti-correlations between the zero wave 
vector and the pattern Fourier modes. Finally only the 
eigenmodes with maximum damping contribute to the 



ffuctuations of the total transverse momentum, there- 
fore the total transverse momentum has the least possible 
ffuctuations. 

Our results are obtained from both semianalytical cal- 
culations based on linearization around the hexagonal 
pattern and from numerical simulations of the nonlin- 
ear system. Some of our results and predictions are very 
general, and depend only on basic symmetry properties 
of the system such as translational invariance. Thus we 
expect that similar structure and properties of the ffuc- 
tuations and correlations can be found in other nonlinear 
systems displaying hexagonal patterns. 

From a computational point of view, correlations in the 
far ffeld intensity peaks are practically independent of the 
system size and therefore can be calculated accurately in 
relatively small systems, provided all the relevant har- 
monics of the pattern are considered. Good statistics 
can be obtained integrating the nonlinear equations over 
relatively short times (even though the near ffeld ffuc- 
tuations are quite far away from reaching a stationary 
value). Of course, alternatively a linear semianalytical 
approach as the one described in Sec. 
used to calculate far-ffeld correlations. 
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